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1.  SUMMARY 


Recent  laboratory  and  theoretical  studies  have  shown  that  the  Green  functions  of  a 
structure  can  be  obtained  from  the  cross-correlation  of  diffuse  wavefields.  The  idea  has 
now  found  rapid  applications  in  seismology,  in  particular,  surface  waves  have  been  found 
to  be  most  easily  retrievable  from  the  cross-correlations  of  seismic  coda  or  ambient 
noise.  Dispersion  measurements  made  on  the  estimated  Green  functions  present 
significant  advantages  over  traditional  measurements  that  are  based  on  earthquakes  with 
limited  distribution.  The  method  is  particularly  useful  in  surface- wave  path  calibration 
and  tomographic  mapping  in  aseismic  regions  and  at  short  periods  (below  20  s).  We 
have  recently  performed  ambient  noise  tomography  of  China  using  the  new  national 
China  Seismic  Network  and  surrounding  global  and  regional  stations.  For  most  of  the 
station  pairs,  we  can  retrieve  very  good  Rayleigh  waves  from  ambient  noise  correlations 
using  18-months  of  continuous  data  at  all  distance  ranges  (100  km  to  5000  km)  and  for 
periods  down  to  about  8  s.  The  combination  of  stations  achieves  good  and  relatively 
uniform  coverage  for  most  parts  of  China.  We  obtain  Rayleigh  wave  group  velocity 
dispersion  measurements  using  a  frequency-time  analysis  method  and  construct  Rayleigh 
wave  group  velocity  maps  for  periods  from  8  s  to  60  s.  Our  preliminary  dispersion  maps 
show  significant  features  that  correlate  with  surface  geology.  Two  major  features  stand 
out  in  particular.  First,  the  major  basins  with  thick  sediment  deposition  correlate  well 
with  slow  group  velocities  at  shorter  periods  (10  to  20  s).  The  major  basins,  including 
Tarim,  Junggar,  Qadaim,  Sichuan,  Bohai,  Songliao,  Southern  North  China,  and 
Jianghan,  are  all  well  delineated  by  slow  velocities.  Second,  variations  in  crustal 
thickness  correlate  with  group  velocities  for  periods  around  30  s.  The  major  trend  of  crust 
thickening  from  the  east  to  west  is  well  represented  by  the  velocity  decreases  from  east 
to  west.  We  propose  several  methods  to  estimate  the  reliability  and  error  of  Rayleigh 
wave  dispersions  derived  from  the  ambient  noise. 
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2.  INTRODUCTION 


The  overall  objective  of  this  project  is  to  obtain  surface  wave  dispersion  measurements 
from  ambient  seismic  noise  correlations  of  the  Chinese  backbone  stations  (CNSN)  and 
global  and  temporary  stations  in  the  surrounding  regions,  and  use  these  measurements  to 
produce  surface  wave  dispersion  maps  of  China.  The  ultimate  goal  is  to  improve 
capability  for  detecting  and  monitoring  small  events  using  the  Ms:  mb  discriminant  by 
improving  path  calibration  of  surface  wave  propagation,  particularly  in  aseismic  areas. 
Below  I  give  a  brief  introduction  to  ambient  noise  tomography.  I  then  describe  the  data 
used  in  this  study  and  the  procedures  of  data  processing  in  Section  3.  Results  are 
presented  in  Section  4  and  concluding  remarks  are  in  Section  5.  This  project  is  in 
collaboration  with  Dr.  Sihua  Zheng  (China  Earthquake  Administration)  and  Dr.  Mike 
Ritzwoller  (University  of  Colorado  at  Boulder).  The  latest  results  have  been  presented  in 
the  AGU  2007  Fall  Meeting  as  an  invited  paper  (Song  et  al.,  2007). 

Recent  theoretical  and  laboratory  studies  have  shown  that  the  Green  functions  of  a 
structure  can  be  obtained  from  the  cross-correlation  of  diffuse  wavefields  (Weaver  and 
Lobkis,  2001;  Lobkis  and  Weaver,  2001;  Derode  et  al.,  2003;  Snieder  2004;  Wapenaar, 
2004;  Larose  et  al.,  2004;  see  also  reviews  by  Larose  et  al,  2006  and  Campillo,  2006). 
The  basic  idea  is  that  linear  waves  preserve,  regardless  of  scattering,  a  residual  coherence 
that  can  be  stacked  and  amplified  to  extract  coherent  information  between  receivers 
(Weaver,  2005).  The  basic  approach  can  be  traced  back  much  earlier  studies  of  random 
fields  in  seismology  (Aki,  1957;  Toksoz,  1964;  Claerbout,  1968),  in  helioseismology 
(Rickett  and  Claerbout,  1995),  and  in  acoustics  (Eagle,  1981;  Weaver,  1982).  The  idea 
has  now  found  rapid  applications  in  seismology.  In  particularly,  surface  waves  have  been 
found  to  be  most  easily  retrievable  from  the  cross-correlations  of  seismic  coda  (Campillo 
and  Paul,  2003;  Paul  et  al.,  2005)  or  ambient  noise  (Shapiro  and  Campillo,  2004;  Shapiro 
et  al.,  2005;  Sabra  et  al.,  2005a,  2005b)  between  two  stations,  although  a  retrieval  of  the  P 
waves  has  been  reported  (Roux  et  al.,  2005). 
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Tomographic  studies  of  surface  waves  from  ambient  noise  correlations  were  first 
conducted  by  Shapiro  et  al.  (2005)  and  Sabra  et  al.  (2005b)  in  S.  California.  The  results 
were  quite  impressive,  showing  nice  correlations  between  surface  wave  speeds  and  major 
geological  units.  The  approach  has  now  been  rapidly  applied  in  different  regions  of  the 
world:  Soutwest  China  (Yao  et  al.,  2006),  Europe  (Yang  et  al.,  2007),  and  S.  Korea 
(Kang  and  Shin,  2006;  Cho  et  al.,  2006;  Yoo  et  al.,  2006). 

The  ambient  noise  tomography  overcomes  some  of  the  important  limitations  of 
conventional  methods  based  on  energetic  sources  (mainly  earthquakes).  (1)  The  data 
coverage  is  limited  by  source  distribution.  Earthquake  distribution  is  very  un-uniform, 
resulting  in  very  un-uniform  sampling.  This  is  particularly  a  problem  in  aseismic  regions. 
(2)  Tomography  requires  knowledge  of  earthquake  location  and  perhaps  focal 
mechanisms.  The  knowledge  is  not  accurate.  (3)  Constrained  by  source  distributions, 
higher  frequency  signals  are  attenuated  by  intrinsic  attenuation  and  the  scattering  of  the 
media  when  the  path  is  long,  limiting  the  resolution  of  seismic  imaging. 

3.  DATA  AND  METHOD 

The  data  we  used  include  the  recent  Chinese  backbone  stations  and  global  stations  in  the 
surrounding  regions  (Fig.  1).  The  Chinese  backbone  network  consists  of  48  stations, 
which  are  distributed  fairly  uniformly  in  China.  We  used  47  of  the  48  stations,  because 
one  station  (Xisha)  has  no  data.  The  network  has  been  referred  to  as  China  National 
Seismic  Network  (CNSN),  China  National  Digital  Seismic  Network  (CNDSN),  Center  of 
China  Digital  Seismic  Network  (CCDSN)  stations,  and  China  Seismic  Network  (CSN). 
We  refer  here  as  China  National  Seismic  Network  (CNSN).  Except  for  stations  of  the 
China  Digital  Seismic  Network  (CDSN),  which  was  established  between  China/State 
Seismological  Bureau  and  the  US  Geological  Survey  in  1986,  most  of  the  new  stations 
started  operating  in  early  2000.  The  bandwiths  of  the  instruments  are  from  20  Hz  to  over 
120  s.  Recently,  25  CNSN  stations  (including  10  CDSN  stations)  have  been  designated 
as  international  exchange  stations.  Data  access  to  the  rest  (23)  of  the  CNSN  stations  is 
more  limited. 
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Fig.  1.  Distribution  of  China  Seismic  Network  stations  and  global  stations  used  in  this 
study. 

We  use  techniques  and  software  that  are  developed  by  Ritzwoller  and  co-workers.  The 
basic  techniques  have  been  discussed  in  great  detail  by  Bensen  et  al.  (2007).  Below  is  an 
outline  of  the  data  processing  and  imaging  steps  we  use.  The  first  step  is  data  pre¬ 
processing,  cross-correlation,  and  stacking.  The  pre-processing  includes  clock 
synchronization  (some  Chinese  stations  use  local  time  instead  of  GMT),  removal  of 
instrument  response,  time-domain  filtering,  temporal  normalization  and  spectral 
whitening.  The  purpose  is  to  reduce  the  influence  of  earthquake  signals  and  instrument 
irregularities  and  to  enhance  the  strength  and  bandwith  of  ambient  noise  in  the  cross¬ 
correlations.  The  result  of  the  first  step  is  the  empirical  Green  function  from  the  cross¬ 
correlation.  Second,  if  the  signal-to-noise  ratio  (SNR)  is  large  enough,  group  speed 
curves  are  measured  in  each  period  band  based  on  a  frequency-time  analysis  (FT AN) 
(Ritzwoller  and  Levshin,  1998).  Third,  repeat  dispersion  measurement  from  ambient 
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seismic  noise  for  different  time  periods  and  use  the  result  to  estimate  measurement 
uncertainties.  Similarly,  if  there  are  dense  arrays  of  stations  as  in  North  China,  S.  Korea, 
and  Japan,  we  can  repeat  dispersion  measurements  between  a  remote  station  and  stations 
that  are  close  by  to  estimate  measurement  uncertainties.  Fourth,  the  final  stage  of  ambient 
noise  tomography  is  group  speed  tomography  itself.  Dispersion  measurements  between 
station-pairs  can  be  used  to  produce  group  speed  maps  exactly  as  earthquake  data  are 
used. 


4.  RESULSTS  AND  DISCUSSION 

4. 1.  Ambient  Noise  Correlations  and  Dispersion  Measurements 

Our  analysis  is  based  on  a  total  of  18  months  of  continuous  data  that  we  have  collected. 
The  data  cover  the  periods  from  Nov.  2003  through  Oct.  2004  and  from  Jan.  2007 
through  June  2007.  For  most  of  the  station  pairs,  we  can  retrieve  very  good  Rayleigh 
waves  from  ambient  noise  correlations.  Fig.  2  shows  typical  examples  of  empirical  Green 
functions  of  Rayleigh  waves  retrieved  from  ambient  noise  correlations.  The  cross¬ 
correlations  show  strong  arrivals  at  different  settings  (near  the  coast  or  well  into  the 
continental  interior)  and  at  both  relatively  low  frequencies  (20-50  s)  and  high  frequencies 
(5-20  s). 
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Fig.  2.  Ambient  noise  cross-correlations  relative  to  two  CNSN  stations:  one  near  the 
coast  (SSE,  Shanghai)  (upper  panels)  and  one  in  the  continental  interior  (GTA,  Gangsu 
Province)  (lower  panels).  The  plot  shows  cross-correlations  of  the  station  with  all  the 
other  CNSN  stations  using  1  year  of  data.  The  cross-correlations  are  filtered  at  two 
frequency  bands  (5  to  20  s  and  20  to  50  s). 


For  pairs  with  Rayleigh  wave  SNR>10,  we  measured  group  velocity  dispersion  curves 
using  the  FTAN  method  (Ritzwoller  and  Levshin,  1998).  The  measurement  is  very  stable. 
Fig.  3  shows  two  examples  of  Rayleigh  wave  group  velocity  dispersion  measurements 
from  the  empirical  Green  function  retrieved.  Note  the  group  velocity  measurements  can 
extend  to  periods  of  10  s  or  shorter  even  for  station  pairs  that  are  separated  over 
thousands  of  kilometers.  The  AXX-QIZ  path  shows  group  velocities  very  similar  to  the 
predictions  from  a  global  3D  model  based  on  earthquake  data  (Shapiro  and  Ritzwoller, 
2002).  However,  the  HTA-BRYK  path,  which  samples  the  Tarim  Basin,  shows 
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significantly  lower  group  velocities  at  short  periods  (below  30  s).  This  is  consistent  with 
thickness  sediments  of  the  basin  (see  discussion  below). 
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Fig.  3.  Examples  of  station- station  cross-correlations  and  frequency-time  analysis 
(FTAN)  (Ritzwoller  and  Levshin,  1998)  to  retrieve  dispersion  curves  (blue  lines)  for 
Rayleigh  waves  of  station  pair  AXX-QIZ  (upper  panels)  and  HTA-BRVK  (lower  panels). 
From  left  to  right  are  map  of  the  ray  path,  the  Green  function  filtered  at  different  periods, 
and  FTAN  plot,  respectively.  In  the  FTAN  plot,  the  black  curve  is  the  prediction  from  the 
3-D  model  by  Shapiro  and  Ritzwoller  (2002),  which  is  used  for  phase-matching  filter  in 
the  data  analysis  and  for  comparison  with  actual  measurement  (in  blue). 


We  have  obtained  thousands  of  dispersion  measurements  with  SNR>10  for  periods  8  s  to 
70  s  (Fig.  4).  The  SNR  is  the  signal-to-noise  ratio  defined  as  the  peak  amplitude  of  the 
Rayleigh  wave  divided  by  the  rms  amplitude  of  the  background  noise.  The  best  observed 
frequency  band  is  10  to  30  s  with  retrievals  of  50  to  80%  of  the  possible  pairs. 


7 


2000- 


C/D 

HH-1 

C 

0 

E 

0 

v_ 

=3 

ca 

o 

E 


0 

_Q 

E 

=j 

c 


total  number  of  paths 


1  500  - 


1000  - 


500  -  * 


■ 

*•••.•  »J 8  month  stacks  {SNR  >10) 

•  • 

*•  •  • 

*  *  *••  I  **  • 

1 2  month  stacks  (SNR  >  1 0)  *  *  -  *  l 


10  20  30  40  50 

period  (s) 


60 


70 


Fig.  4.  Number  of  dispersion  measurements  obtained  for  different  periods  in  this  study. 
The  total  number  of  possible  paths  (station  pairs)  is  about  1700.  The  numbers  of 
measurements  obtained  from  18-months  of  continuous  data  are  in  red  dots.  The  numbers 
show  about  20%  of  improvement  over  the  numbers  of  measurements  obtained  using  only 
12-months  of  data  (blue  dots).  Rayleigh  waves  are  retrieved  over  a  wide  frequency  band 
(8  to  70  s).  The  best  observed  frequencies  are  10  to  30  s. 


4.2.  Surface-Wave  Tomography 


Using  the  dispersion  measurements  we  have  made,  we  obtained  tomographic  maps  of 
group  velocities  at  frequencies  from  8  s  to  60  s  (Fig.  5).  The  results  show  some 
remarkable  features  that  correlate  with  large-scale  geological  structures  of  China.  Major 
basins  are  well  delineated  with  low  velocities  at  periods  8  s  and  20  s,  including  Bohai 
Basin  (North  China  Basin),  Sichuan  Basin,  Qaidam  Basin,  and  Tarim  Basin.  The  Yangtz 
Craton  also  shows  up  well  with  high  velocities.  At  a  longer  period  (25-50  s),  the  group 
velocity  maps  display  striking  bi-modal  distribution  with  high  velocity  in  the  east  and 
low  velocity  in  the  west,  which  corresponds  very  well  the  thinner  crust  in  the  east  and 
much  thicker  crust  in  the  west  (e.g.  Liang  et  al.,  2004).  The  NNE-SSW  trending 
boundary  between  fast  and  slow  velocities  (around  longitude  108°E)  coincides  with  the 
sharp  topography  change  and  with  the  well-known  Gravity  Lineation. 

Fig.  6  shows  side-by-side  comparison  between  the  group  velocity  map  at  15  s  and 
sediment  thickness.  Fig.  7  shows  comparison  between  the  group  velocity  map  at  30  s  and 
the  crustal  thickness.  The  correlations  are  quite  remarkable  if  we  compare  the  group 
velocities  along  a  certain  profile  of  interest  (Fig.  8). 
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Fig.  5.  Rayleigh  wave  group  velocities  obtained  in  this  study.  Show  are  the  maps  at 
periods  8,  15,  20,  25,  30,  40,  50,  60  s,  respectively. 
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Fig.  6.  Comaprison  of  group  velocity  map  at  15  s  (left)  and  sediment  thickness  (right). 
Plotted  in  the  background  are  major  block  boundaries  and  basin  outlines  (Liang,  Song 
and  Huang,  JGR,  2004).  The  major  basins  (including  Tarim,  Junggar,  Qadaim,  Sichuan, 
Bohai,  Songliang,  Southern  North  China,  Jianghang)  are  well  delineated  by  slow 
velocities. 
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Fig.  7.  Comparison  of  group  velocity  map  at  30  s  (left)  and  crustal  thickness  (right).  The 
major  trend  of  crust  thickening  from  the  east  to  west  (right)  is  well  represented  by  the 
velocity  decreases  from  east  to  west  (left). 
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A) 


profile  from  80/38  to  125/38 


profile  from  80/38  to  125/38 


Fig.  8.  Comparison  of  group  velocities  and  sediment  thickness  and  crustal  thickness 
along  the  selected  path.  (A)  The  selected  profde  along  latitude  38  deg  passes  three  major 
basins  (Tarim,  Ordos,  and  Bohai).  (B)  The  three  major  basins  correlate  well  with  slow 
velocities  at  periods  10  to  20  s.  The  surface  elevation  (middle,  shaded  green)  is  plotted 
for  reference.  The  elevation  and  sediment  thickness  are  plotted  using  the  same  scale.  (C) 
The  general  trend  of  decreasing  crustal  thickness  from  west  to  east  is  well  represented  by 
increasing  group  velocities  at  30  s.  However,  the  group  velocity  shows  more  structure 
than  the  smooth  crustal  thickness  curve  from  the  global  reference  model  (CRUST  2.0), 
suggesting  likely  more  complex  Moho  structure.  At  period  50  s,  the  trend  is  not 
observable  any  more  as  the  surface  waves  sample  deeper  into  the  mantle. 
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4.3.  Robustness  of  Measurements 


For  a  diffuse  field,  it  has  been  proved  that  cross-correlation  yields  the  Green  function 
between  the  stations,  as  discussed  in  the  introduction.  If  the  field  is  not  diffuse,  how  can 
we  validate  that  the  Rayleigh  waves  we  have  obtained  from  the  cross-correlation  are  the 
true  solution  of  the  media.  If  it  is  not  a  true  solution,  what  is  the  error  of  the  estimate? 
Some  ideas  have  been  proposed  before.  (1)  Comparing  the  cross-correlation  with  surface 
wave  generated  by  an  earthquake  along  the  same  path.  (2)  Comparing  correlations  using 
data  from  different  seasons  of  the  year.  Because  there  will  be  different  noise  sources  from 
the  oceans,  the  consistency  of  the  correlations  gives  a  measure  of  stability  and  error  of  the 
Green  function. 

We  propose  further  validation  methods.  (3)  Model-based  validation.  The  surface  wave 
tomographic  maps  can  be  used  to  test  how  good  the  models  are.  In  this  regard,  the 
discussion  above  on  comparing  the  tomographic  maps  with  geological  structure  and 
processes  is  particularly  relevant.  A  good  model  should  be  consistent  with  well-know 
geological  features.  (4)  Temporal  comparison:  Comparing  correlations  with  different 
sliding  time  windows  (Fig.  9).  For  noise  sources  that  may  not  be  seasonal,  the  method 
provides  a  way  to  measure  the  variations  of  different  time  windows.  The  dispersion 
curves  are  very  stable.  Longer  time  windows  improve  the  stability  at  longer  periods,  but 
also  at  the  shorter  periods.  (5)  Spatial  comparison:  Comparing  correlations  for  pairs  along 
similar  paths  (Fig.  10).  Green  functions  between  a  station  to  two  close  stations  should  be 
similar  as  they  sample  similar  structure.  This  is  indeed  the  case  for  YCH  to  GZH  and 
YCH  to  SZN.  YCH  falls  on  the  great  circle  path  passing  GZH  and  SZN,  which  are 
relatively  close  to  each  other. 

5.  CONCLUSIONS 

We  have  obtained  Rayleigh  wave  group  velocity  measurements  using  ambient  noise 
correlations  of  18  months  of  continuous  data  of  CNSN  and  some  global  stations.  The 
Rayleigh  waves  can  be  retrieved  from  about  8  s  to  70  s.  The  best  observed  frequency 
band  is  10  to  30  s  with  retrieval  rate  of  50  to  80%  of  the  pairs.  Rayleigh  group  velocity 
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tomographic  maps  show  remarkable  correlations  with  major  tectonic  features  of  china. 
The  correlations  between  the  tomographic  models  and  the  tectonic  features  provide  an 
important  validation  of  the  methodology,  the  dispersion  measurements,  and  the 
tomographic  results.  Studies  of  the  temporal  and  spatial  consistency  also  suggest  that  our 
measurements  are  robust.  Further  studies  will  quantify  the  errors  of  the  measurements 
and  tomographic  models. 


12-monlh  slacks(l9),  step  10  days 


6-month  stacks(l3),  step  1  month 


Dispersion  of  Group  Velocity 

(GTA_OIZ) 


Dispersion  of  Group  Velocity 

(GTA_TATO) 


Fig.  9.  Comparison  of  dispersion  curves  from  different  time  windows.  Two  examples  are 
BJT  and  BRVK  pair  (upper  panels)  and  GTA  and  QIZ  pair  (lower  panels).  For  each  path 
from  left  to  right  are  the  map  of  the  path,  dispersion  curves  obtained  from  12-months  of 
data  with  a  sliding  time  window  of  10  days  (total  of  19  curves),  and  dispersion  curves 
obtained  from  6-months  of  data  with  a  sliding  time  window  of  1  month  (total  of  13 
curves). 
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Dispersion  of  Group  Velocity 

(YCH_GZH{R)/SZN(B)) 


period  (s) 


Fig.  10.  Comparison  of  dispersion  curves  between  YCH  and  GZH  (red)  and  between 
YCH  and  SZN  (blue).  The  two  paths  are  very  similar.  The  group  velocities  are  also  very 
similar  for  periods  where  they  overlap  (10  to  40  s). 
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